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We briefly review Boltzmann-Gibbs and nonextensive statistical mechanics as well as their con- 
nections with Fokker-Planck equations and with existing central limit theorems. We then provide 
some hints that might pave the road to the proof of a new central limit theorem, which would play 
a fundamental role in the foundations and ubiquity of nonextensive statistical mechanics. The basic 
novelty introduced within this conjectural theorem is the generalization of the hypothesis of indepen- 
dence of the random variables being summed. In addition to this, we also advance some nonlinear 
dynamical (possibly exact) relations which generalize the concepts of Lyapunov exponents, entropy 
production per unit time, and their interconnection as first proved by Pesin for chaotic systems. 



I. INTRODUCTION 

As well known, thermodynamics is the basic branch 
of physics which focuses on the generic connections be- 
tween variables (temperature, pressure, volume, energy, 
entropy and many others) that play an important role 
in the description of the macroscopic world. Boltzmann 
and Gibbs provided a magnificent connection of thermo- 
dynamics with the microscopic world lj ,2j. This connec- 
tion, normally referred to as Boltzmann-Gibbs (BG) sta- 
tistical mechanics (or simply statistical mechanics since 
it was basically the only one to be formulated along more 
than one century), turns out to be the appropriate one 
for ubiquitous systems in nature. It is based on the fol- 
lowing axiomatic expression for the entropy: 
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where pi is the probability associated with the i micro- 
scopic state of the system, and k is Boltzmann constant. 
In the particular case of equiprobability, i.e., pi = 1/W 
(V«), Eq. (1) yields the celebrated Boltzmann principle 
(as referred to by Einstein himself 0): 



Sbg = klnW . 



(3) 



For continuous variables, the BG entropy is written as 
Sbg^ - J dxp{x) In p{x) (xeM'^), (4) 



with 



dxp{x) — 1 



(5) 



If X happens to carry physical units, we write Eq. (4) as 
follows: 



Sbg = 



dx p{x) \ny{^^ar^p{x) 



(6) 



where (7^ > carries the same units as the r"* component 
of the d-dimensional variable x. Clearly, when x carries 
no units (i.e., when x G E**), we take ct^ = 1 (Vr). In 
fact, everytime this is possible (and it is possible most of 
the time), we shall adapt the physical units in such a way 
that CTr = 1 (Vr) even when x does have physical units. 
Consistently, unless otherwise specified, we shall use Eq. 
(4) for the general continuous case. If we consider the 

particular case p{x) = X^i^i Pi^i^ — ^i) (with Pi — 
1, {xi} being some set of values, and d{z) being Dirac' s 
delta distribution), Eq. (4) recovers Eq. (1). 

For quantum systems, the BG entropic form is written 
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with 



From now on, and without loss of generality, we shall 
take k equal to unity. 
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Trp = 1 , 



(8) 



p being the density operator or matrix. When the W x 
W matrix p is diagonalized, it shows the set {pi} in its 
diagonal. In what follows, depending on the context, we 
shall use either the discrete form (Eqs. (1) and (2)), or 
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the continuous form (Eqs. (4) and (5)), or the matricial 
form (Eqs. (7) and (8)). 

In spite of its tremendous power and usefulness, the 
BG concepts and statistical mechanics appear to be not 
universally applicable. Indeed, there is a plethora of nat- 
ural and artificial systems (see, for instance, 4] and ref- 
erences therein) for which they do not provide the ad- 
equate mathematical frame for handling physically rel- 
evant quantities. This fact started being explicitly rec- 
ognized at least as early as in 1902 by Gibbs himself: 
see page 35 of 0, where he addresses anomalies related 
to systems such as gravitation. A formalism becomes 
therefore desirable which would address such anomalous 
systems. A vast class of them (although surely not all of 
them) appears to be adequately discussed within a gen- 
eralization of the BG theory, frequently referred to as 
nonextensive statistical mechanics. This theory was first 
introduced in 1988 @, and then refined in 1991 '6^ and 
1998 0. It is based on the following generalization of 
Sbg'- 

S, = ^^=^^' (q e K; 5i = Sbg) ■ (9) 

Expressions (4), (6) and (7) are respectively general- 
ized into 



^ 1 - J dx [p{x)]'i 



(10) 



Levy-Gnedenko central limit theorems, and argue about 
the possible formulation of a new theorem that would 
generalize the standard limit theorem. Such a theorem 
would provide an important mathematical cornerstone 
for nonextensive statistical mechanics and its ubiquity in 
nature. Providing hints that could help formulating and 
proving the theorem constitutes the main reason of the 
present paper. 



II. STATISTICAL MECHANICS 

Entropy is necessary to formulate statistical mechan- 
ics but it is not sufhcient. Indeed, we must also intro- 
duce the concept of energy. The easiest (and more fre- 
quently used) way to do so is addressing the so called 
(by Gibbs) canonical ensemble. It corresponds to the 
ubiquitous physical situation in which the system of in- 
terest is in contact with a (large by definition) thermo- 
stat which, at equilibrium (or at the physically relevant 
stationary state, more generally speaking), imposes to 
the system its temperature. The system is typically de- 
scribed by a quantum Hamiltonian, and is characterized 
by the spectrum of energies {Ei} (i = 1, 2, W) defined 
as the eigenvalues associated with the Hamiltonian and 
its boundary conditions. The probability distribution at 
the relevant stationary state is the one which extremizes 
the entropy under the norm and the energy constraints. 
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For equiprobability (i.e., pi = 1/W, Vi), Eq. (9) yields 
Sq = lnqW , (13) 
with the q-logarithm function defined as 

z^-'i - 1 



A. Boltzmann-Gibbs statistical mechanics 

At thermal equilibrium we must optimize Sbg 
given by Eq. (1)) with the norm constraint as given by 
Eq. (2), and the energy constraint given as follows: 
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'^PiEi = Ubg , 



(16) 



where U bg is given and referred to as the internal en- 
ergy. By following the Lagrange method, we define the 
quantity 



In^z : 
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(z e M; z > 0; Iniz = Inz) . (14) -^bg ^ SBG + a{^p^-lj - ^3{^p^E,-UBG^ , (17) 
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Its inverse function, the q- exponential, is given by 

e^^[l + (l-g)z]i/(i-^) {el = e^) (15) 

if the argument 1 -I- (1 — q)z is positive, and equals zero 
otherwise. 

In Section II we briefly review a few known results 
concerning the probability distributions which extremize 
the entropy, and are ultimately associated with macro- 
scopically stationary states. In Section III we briefly 
review available results related to Fokker-Planck equa- 
tions. In Section IV we review the standard and the 



where a and (3 are the Lagrange parameters (their signs 
have been chosen following tradition). The extremizing 
condition S^bg/^Pj — yields 

=e"-i"^^^ {j = l,2,...,W). (18) 

The use of Eq. (2) allows the elimination of the param- 
eter a. We then obtain the celebrated Boltzmann-Gibbs 
weight 



Zbg 



(z = l,2,...,I^), 



(19) 
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where (3 is connected with the thermostat temperature T 
through [3 = l/T, and the partition function is defined 
as follows: 



w 



(20) 



It is precisely the present Eqs. (19) and (20) that are 
referred to as dogma by the mathematician Takens! HI . 

Clearly, if we replace the probability distribution (19) 
into Eq. (16), we obtain the thermodynamically impor- 
tant relation between the inverse temperature f3 and the 
internal energy Ubg- 

As a subsidiary comment, whose relevance will become 
transparent later on, let us remark that the BG weight 
(19) can be seen as the solution of the linear ordinary 
differential equation 



dy 

dx 



ay (2/(0) - 1) 



Indeed, its solution is given by 



(21) 



(22) 



which reproduces Eq.(19) through the identification 
{x,a,y) = {Ei,-(3,ZBGPt)- 



B. Nonextensive statistical mechanics 

We want now to optimize Sq (as given by Eq. (9)) 
with the norm constraint still given by Eq. (2), and the 

I): 



energy constraint generalized as follows (see 



Z^i=l Pi 



(23) 



(referred to as q-expectation value or q-mean value) or, 
equivalently, 



W 



(24) 



The functional (17) is generalized into 
w w 



Sq + a^p,-l]-P^pl{E,-Uq)\, (25) 

i=l i=l 

and the extremizing condition S^q/Spj = yields 

= {J = 1,2,...,W). (26) 

The use of Eq. (2) allows, as before, the elimination of 
the parameter a, thus obtaining the generalized weight 



-t3{Ei-Ug) 



ii^l,2,...,W), (27) 



with 
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(28) 



This probability distribution corresponds to a maximum 
[minimum) of Sq for g > (q < 0). For g = 0, the en- 
tropy is constant, namely 5*0 = — 1, and the distribu- 
tionisgivenbyK = [l~(i{E^~Uo)y Y.J=i[^-P{E,~U^)] 
(we recall the cutoff of the q-exponential function for 
q < 1, i.e., the states for which 1 — P{Ei — Uq) < 
do not contribute). 

As a subsidiary comment, let us remark that the weight 
(27) can be seen as the solution of the nonlinear ordinary 
differential equation 



dy_ 

dx 



Indeed, the solution is given by 
2/ = er , 



(29) 



(30) 



which reproduces Eq.(27) through the identification 

{x,a,y) = {Ei -Uq,-(3,ZqPi). 



C. The Gaussian case 

Let us now assume a stochastic real variable x such 
that 



(x) 



and 



dxxp{x) — , 



/oo 
dx x'^p(x) = D . 
-oc 



(31) 



(32) 



where D > is given. The maximization of entropy 
(4) with constraints (5) and (32) yields the celebrated 
Gaussian distribution, more precisely 



p{x) = \ — e 

V TT 



(33) 



where /? is a Lagrange parameter. By using Eq. (32), we 
immediately obtain /3 ~ 1/21?, hence 



p{x) 



y/2nD 



(34) 



We straightforwardly verify that constraint (31) is satis- 
fied as well. 

A more general case is to assume the following con- 
straints: 



dxxp{x) — C . 



(35) 
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and 

/oo 
dx {x- {x)fp{x) ^ D . 

The extremizing distribution is then given by 
p{x) = 

or, equivalently, 



p{x) 



-{x-Cr/{2D) 



(36) 



(37) 



(38) 



D. The g-generalization of the Gaussian case 

Let us assume that the constraints are now even more 
general, namely 



and 



^ f^^dxx[p{x)]i _ 
C^dx[p{x)]^ 



_ jr^dxix~{x),nPi^)]^ 



C^dxipix)]" 



(39) 



Dg, (40) 



where P{x) = [p{x)]'^ / J dy[p{y)]'^ is, in the literature, 
referred to as the escort distribution. 

The associated distribution extremizing Sq is now 
given by 



p{x)^Aq^e;^^^-^^^^^' 



(41) 



where we have flQ] A, = ^y{q - l)/7rr(l/(q-l))/r((3- 
q)/[2/iq ~ 1)]) for q > 1, and - ^{1- <z)Ar((5 - 
3q)/[2(l - q)])/r{{2 - q)/{l - q)) for q < 1, r(^) being 
the Riemann function. 

The use of constraint (40) straightforwardly provides 
/3 = l/[(3 — q)Dq], which, replaced in Eq. (41), yields 

(42) 

In the q < 1 case, the support is compact, and the 
distribution vanishes outside the interval |a; — Cq| < 
^{3-q)Dq/{l-q) 
If (x 



q — 0, distributions (41) and (42) take respec- 
simple forms 

(43) 



tively the simple forms 

p{x) ^ Aqy^e-/^' 
and 

A, 



p{x) 



y/{3~q)Dq, 



^xy[{3-q)D,] ' 



(44) 



which appear frequently in various contexts. 

It can be seen ^^7] that these distributions are ana- 
lytic extensions (to real values of q) of the Student's t- 
distribution and the r-distribution, for q > 1 and q < 1 
respectively. Consistently, there is an asymptotic long 
power-law tail for q > 1, and a compact support for q < 1. 
There is an upper bound for q, namely q = 3, imposed by 
the norm constraint (5). More precisely, the admissible 
values of q are q < 3. It deserves to be mentioned that 
also constraint (40) has an upper bound in order to be 
finite, which happens to be precisely the same, i.e. q = 3. 



DIFFUSION AND FOKKER-PLANCK 
EQUATIONS 



III. 



Normal diffusion (i.e., the one which satisfies (x^) cx i, 
typical of Brownian motion 12]) is characterized by the 
heat equation (the simplest form of the Fokker-Planck 
equation) 



dp{x,t) j^d'^p{x,t) 



{D>0). 



dt dx^ 
By assuming, at t = 0, the paradigmatic form 

p{x,0) = S{x) , 

we obtain the following exact solution: 



p{x,t) 



V27rL>t 



(t>0) 



(45) 



(46) 



(47) 



There is of course an infinity of possible generaliza- 
tions of Eq. (45), which is linear and defined through 
integer derivatives. We address here two important such 
generalizations, both of them associated with anoma- 
lous diffusion, i.e., violating the relation {x"^) cx t The 
first one remains linear but replaces the second deriva- 
tive of the right term by a fractional derivative; it yields 
(x^) — > oo (yt > 0). The second one is nonlinear but pre- 
serves the integer derivatives; it yields (x^) cx (a 7^ 1). 

The first one is as follows: 

dpix,t) _ ^dMxJl (^>o.o<^<2) (48) 



dt 



dx'' 



The exact solution associated with the t — condition 
(46) is given by 



p{x,t) 



1 



—L,{xl{Dtf''<) (t>0), (49) 



where L^{z) is the Levy distribution of index 7. The 
7 — > 2 limit obviously corresponds to the Gaussian solu- 
tion (47). The 7 = 1 particular case corresponds to the 
Cauchy-Lorentz distribution 



. X , 1 Dt 

p{x,t) = Li{x/{Dt)) = -- 



TT (Dt)^ + x2 



(50) 
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For all other values of 7 different from 1 and 2, the Levy 
distribution has no direct analytic expression, and is ex- 
pressed only through its Fourier transform. 

The second generalization of Eq. (45) is as follows: 



dt 



The exact solution associated with the t 
(46) is given by 



). (51) 
condition 



with 



-x^/[{J.-q)(Dtf'^''-'i'>] 



{t>0) 
(52) 



q = 2-iy <3, 

with D > if > 0, and D < if < 0. 
the diffusion coefficient as D = D/iy, Eq. 
rewritten as follows: 



(53) 

If we rewrite 
(51) can be 



D 



dx^ 
0, becomes 



D- 



[liip{x,t)] , 



(54) 



(55) 



dpjx, t) ^ 
dt 

which, in the limit u - 

dp{x,t) 
dt 

whose solution is known to be the Cauchy-Lorentz dis- 
tribution. 

If we start from any t — Q distribution different from 
(46), the solution p(a;, t) does not coincide with Eq. (52), 
but, nevertheless, it asymptotically approaches IT3 Eq. 
(52) for t 00. In other words, the solution (52) consti- 
tutes an attractor in the space of the distributions, i.e., it 
is a solution which is robust. For all q < 3, x scales [lj| 
with t^/^^"*?). Consequently, in all cases for which the 
second moment of p{x,t) is finite (i.e., for q < 5/3) we 
obtain (a;^) cx ^2/(3-9) ^^Yla.t happens for < 5/3 < 3 is 
that the prefactor of diverges. In this case, it is 

convenient to focus on the width of the distribution, or, 
equivalently, on {x'^)q. 

To close these remarks about Fokker-Planck equations, 
let us write down a quite general one, namely 



d'p{x,t) ^ d''[p{x,t )] 
d\t\^ d\x\ 



2-q 



((<5,7,g)GM'), (56) 



where we have used, for convenience, q rather than v 
(related through Eq. (53)). The regions we address are 
0<(5<1, 0<7<2 and q < 3. Not surprising, 
the solutions for arbitrary {S,^,q) are not known. But, 
nevertheless, we depict the 6=1 regions in Fig. 1. This 
will clarify the particular cases that are related to the 
central limit theorems we are interested in. 



IV. CENTRAL LIMIT THEOREMS 

In one of its simplest versions, the standard central 
limit theorem may be formulated as follows. Let {xi} 



r 

2 




conjectured CLT 






"^^Gaussian CLT 


1- 




'^^Levy-Gnedenko CLT 



2 q 3 



FIG. 1: (5 = 1 regions that are relevant for the Central Limit 
Theorems (CLT). The solution of Eq. (56) for (7,5) = (2, 1) 
is a Gaussian distribution (which optimizes the entropy Sbg, 
and whose variance is finite). For q = 1 and < 7 < 2, the 
solutions of Eq. (56) are Levy distributions (which asymptot- 
ically decay as l/|x|^+^, and whose variance diverges). For 
7 = 2 and q < 3, the solutions are g-Gaussians (which opti- 
mize the entropy Sq, and whose variance is finite for g < 5/3 
and diverges for 5/3 < g < 3; their g- variance is finite for all 
g < 3; they have a compact support for g < 1 and asymp- 
totically decay as l/lxl^/'^"^' for 1 < g < 3). Within the 
stardard CLT, we sum independent variables whose variance 
is finite; within the Levy-Gnedenko CLT, we sum indepen- 
dent variables whose variance diverges; within the conjectural 
CLT, we sum specially correlated variables whose g-variance 



(I = 1, 2, TV) be a set of independent random variables, 
each of them satisfying the same symmetric (with regard 
to X = 0) distribution p{x). Let also p{x) be such that 
(x^) = J^^dx x'^p{x) is finite. Let us define now the 
sum variable 



Z 



N 

1=1 



(57) 



The question to which the central limit theorem answers 
is what is the probability distribution of the random vari- 
able Z when N —^ 00? This answer happens to be 
very simple, namely a Gaussian (in the properly rescaled 
variables). If we call p{Z,N) this distribution (with 
p{Z, 1) = p{x)), to exhibit the Gaussian we must use as 
abscissa Z/^/N^ and as ordinate \/Np{Z,N). Then, we 
gradually see emerging the Gaussian as 00. And, 

what specific Gaussian? The one whose second moment 
precisely coincides with (a;^). 

This theorem clearly is closely related to normal diffu- 
sion, hence to Eq. (45) and to its solution (47), where t 
plays the role of N . 

Let us now address the other central limit theorem 
which is known, namely the Levy-Gnedenko one. In its 
simplest version, it can be formulated as follows. As 
before, the set {xi} (l = l,2,...,iV) is constituted of 
independent random variables, each of them satisfying 
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the same symmetric (with regard to a; = 0) distribu- 
tion p{x). But now, we have the case where (x^) di- 
verges. The question remains the same, namely what 
is the probability distribution of the random variable Z 
when N oo? The answer now is: a Levy distribu- 
tion Ly{Z). With what value for its index 7? All Levy 
distributions decay, for |a;| — > 00, as l/|a;|-'^+'^. This be- 
havior is preserved for the distribution of Z while N 
increases. Let us be more specific: the typical case is 
that for which p{x) decays as with 1 < fj, < 3. 

Then 7 = /i — 1. And what specific Levy distribution 
L-y ? The one which has precisely the same coefficient of 
the asymptotically dominant term. In other words, let us 
consider the finite value lim|2.|^oo[?'(a^)|2;|^]. The particu- 
lar Levy distribution which is asymptotically approached 
when A'^ — !■ 00 is the one which has the same value for 
lim|2.|^oo[^7(2;)|a;|^^'']. To see, on a graphic represen- 
tation, the gradual emergence of the Levy distribution 
while A^ increases, the abscissa must be Z/N^^'^ , and the 
ordinate must be N^^'^p{Z, N). 

This theorem is closely related to the anomalous diffu- 
sion characterized by Eq. (48) and by its solution (49), 
where, as before, t plays the role of A^. 

Let us finally address the conjectural theorem that we 
are focusing in this paper. It is of course the one to be 
associated with the nonlinear Eq. (51) and its solution 
(52). What hypothesis is to be violated in the two preced- 
ing and well studied theorems? The hypothesis of inde- 
pendency! Indeed, we believe that the variables {xi} are 
to be assumed somehow collectively correlated in such a 
persistent manner that the correlation does not disappear 
even in the N 00 limit. For the standard central limit 
theorem, the quantity which is preserved is the second 
moment {x'^) = J_^dx x'^p{x). For the Levy-Gnedenko 
theorem, the quantity which is preserved is the coeffi- 
cient lim|;j,|_,oo[p(a;)|x|^]. For this conjectural theorem, 
some quantity is expected to be preserved. Could it be 

{x^), ^ j:^^ dxxM^)]y jr^ dx [p{x)]'^ ? 

It should be transparently clear at this point that we 
have no definitive arguments for proving this conjectural 
theorem. Various converging paths are nevertheless avail- 
able that might inspire a (professional or amateur) math- 
ematician the way to prove it. Galileo used to say that 
knowing a result is not neglectable in order to prove it! 
It is our best hope that his saying does apply in the 
present case! So, what are these converging paths? Al- 
though naturally intertwingled, let us expose them along 
six different lines. 



1. The q-product hint 

The following generalization of the product operation 
has been recently introduced |[16|. It is called q-product 
and is defined through 

X(»qY =[X^-'^ -^Y^-" -1]^^^^-"^ (qeM.) (58) 



We shall address the case where X > 1 and Y > 1 [T^ . 
This product has the following properties: 

(i) X(E)iY = XY; 

(ii) \ng{X Y) ^lngX + In, Y (whereas \nq{XY) = 
In, A + In, y + (1 - g) (In, X) (In, Y); 

(iii) l/(A®,r)-(l/X)®2-9(l/r); 



(iv) x®,(r®,z) 

{X^-1 + Y^-'J + Z^-i 
1 = A; 



- 2)1/(1-9); 



(v) A(5 

(vi) X(g>gY = Y(g,gX; 

(vii) For fixed {X,Y), X (E)qY monotonically increases 
with q . 

Property (ii) is particularly important since it is di- 
rectly relevant to the extensivity of 5*, that will be ad- 
dressed soon. 

Notice that property (iii) involves, like the solution of 
Eq. (51) with Eq. (53), a q ^ (2 — g) transformation. 
We may apply this product to the number 



W 



eff 



of allowed states in a composed sys- 
tem whose subsystems Ai, A2, Apf have respectively 
Wai,Wa2, ---jWan possible states (by allowed we mean 
that their probability is essentially nonzero). If q = 1 

we have the total number Wai+A2+...+An — Y[i=i 
of states that are not only possible a priori, but even 
generically allowed. But, if g 7^ 1, say g < 1, we expect 
correlations to inhibit (even forbidden occasionaly) some 
of the states, i.e., to be associated with a probability close 
(in some sense) to zero. In this case, the effective total 
eff 

number W^^^^^j^^j^ of allowed states is expected to 

be smaller ih&t Wai+A2+...+An- We expect to have ba- 
sically 



W 



eff 



Ai+A2 + ...+A^ 



Wa, 



N 



1=1 



qWA2 
(1-9) 



>, ... «), Waj. 
-(N-l) 



1/(1-9) 



< WA^+A 



2 + . 



N 

.+a.^IIWa 

1=1 



In particular, for q 
reff 



0, we have 



W 



Ai+A2 + ...+A^ 



N 

1=1 



Wa,]-N + 1. 



(59) 



(60) 



If the A^ subsystems are all equal, we have that 
lnqW^-ff{N) - A In, 11/(1), 



(61) 



the changement of notation clearly being W^ff{N) 



W 



eff 



and W{1) = Wai- Consequently, for 



q = 1, we have 

W^ff{N) = W{N) 
whereas, for q < 1, we have 



[Wil) 



N 



W^ff^N) - |a^ [W{1) 



11-9 



- 1 



1/(1-9) 



< W{N) = [W{1)] 



N 



(62) 



(63) 
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We therefore see that, for g = 1, the number W'^ffiN) of 
nonzero-probabihty states grows exponentially with N, 
whereas, for q < 1, W^ff{N) grows hke a power-law with 
iV, more precisely hke A^i/(i~9). For q = 0, it grows 
linearly with N, more precisely W^ff{N) = N[W{1) — 
1] + 1. 

2. The hint of the extensivity of Sq 

We have recently argued (see and references 
therein) that special correlations may exist between the 
N subsystems Ai,A2,...,Am of a composed system such 
that one (and only one) value of the index q exists which 
ensures the additivity of Sq. The trivial case of course is 
that of Sbg- We consider the case of independency, i.e., 
the joint probabilities given by 



N 



Ai+A2 + ...+A„ 



= X{pf: ,V(zi,Z2,...,z^). (64) 



1=1 



We straightforwardly verify that 



N 



Sbg{Ai +A2 + ... + An) = Y^ Sbg{Ai) 



1=1 



where 



Sbg{Ai+A2 + ... + An) 



(65) 



(66) 



ii=l 12 = 1 ijv = l 



-A2 + ...+AP 



Ai+A2 + ...+A^ 



and 



Wa, 



il=l 

We can also verify 



(67) 



N 



ln[l+(l-q)5g(Ai+A2+...+Ajv)] = J2 ln[l+(l-<z)5g(A0] 

(68) 



1=1 



If the subsystems are all equal, we have 

SBG{N)^NSBGil), (69) 

the notation being self-explanatory. 

If we consider iV = 2 in Eq. (68) we obtain 

Sg{Ai + A2) = Sq{Ai) + Sg{A2) + (1 - q)Sq{Ai)Sg{A2) . 

(70) 

Therefore, Sbg is said to be extensive (or additive). Con- 
sistently, Sq is, unless <? = 1, nonextensive (or nonaddi- 
tive). It is in fact from this property that the statistical 
mechanics we are talking about has been named nonex- 
tensive. As we shall exhibit in what follows, this early 



denomination might (unfortunately) be somewhat mis- 
leading. Indeed, Sq is, for g 1, nonextensive if the 
subsystems are (explicitly or tacitly) assumed indepen- 
dent. But, if they are specially correlated, Sq can in fact 
be extensive for some special value oi q ^ 1 . 

Consider now the following set of joint probabilities 
(clearly corresponding to nonindependent subsystems) : 



A1+A2+.. 
P11...1 



N 



Ai 
Pi 



1=1 



N + 1 



^Ai+A2 + ...+A^ 
Pl...li...l 



(71) 



(72) 



and zero otherwise. Consequently, we have only 
W^lf+^,+...+^„ = (Ell Wa,) -N+1 joint probabih- 
ties which are generically nonzero, consistently with Eq. 
(60). We straightforwardly verify 



N 



So{Ai + A2 + ... + An) = J2 ^o{Ai 



where 

So{Ai+A2 
and 



An) = W 



Ai+A2 + ...-\-A^ 



(73) 



1 , (74) 



SoiAi)=WA,-l (VO 



(75) 



Therefore, for the particular correlations involved in Eqs. 
(71) and (72), it is for g = (and not for q = 1) that Sq 
is additive. 

Unfortunately, we do not know the general explicit 
form of the set of joint probabilities for an arbitrary num- 
ber N of (not necessarily equal) subsystems, correspond- 
ing to q ^ 0, 1. We do know however, for < q < 1, some 
special cases, such as the iV = 2, 3 generic binary subsys- 
tems. They are presented and analyzed in We also 
know, for binary equal subsystems, exact recurrence re- 
lations that enable the calculation of the entire set of 2^ 
joint probabilities associated with N subsystems given 
the entire set of 2^~^ joint probabilities associated with 
A^ — 1 subsystems. This calculation is rather lengthy and 
will be the object of a separate paper. 

Now we present here (in a form which is slightly differ- 
ent jand possibly more transparent, than the one used 
m [18|) the case of TV = 2 equal binary subsystems. 
Consider the set of joint probabilities indicated in Ta- 
ble I. Let us impose the additivity of the entropy, i.e., 
Sq{2) = 25,(1), where 



5,(2) 
and 



l-[p^ + k]" - 2[p (l-p)- K,]" - [(1 - p)2 + k]" 



s,ii) 



q-i 



i-p'^-(i-p)^ 
q-i 



(76) 



(77) 
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The relation between the correlation k, p and q is then 
given by 

2[p« + (l-p)9]-[p2 + K]«-2[j5(l-p)-K]«-[(l-p)2 + K]9 = 1 

(78) 

(With the notation fq{p) = + this relation coin- 
cides with the one indicated in In Fig. 2, typical 
K versus p curves are presented. The lower curve corre- 



sponds to p 



and to (1 — p) + k = 0. The upper 



curve corresponds to p{l — p) — n = 0. The lower curve 
undoubtedly corresponds to q — 0. The upper curve 
could in principle correspond to both cases q — > — oo and 
q oo. Indeed, for n = p{l—p), we have that p'^ + n = p 
and (1 — p)^ + K = 1 — p, hence Sq{2) — Sq{l). Since we 
must also satisfy Sq{2) = 2Sq{l), only two possibilities 
emerge a priori, and these are Sq{l) — (which corre- 
sponds to q ~f oo), and Sq{l) — > oo (which corresponds 
to q ^ — oo). Monotonicity with regard to q suggests 
that the upper curve should correspond to g ^ oo. But 
on the other hand, having, among the four joint prob- 
abilities that are a priori possible in the present N = 2 
system, two zeros (which definitively is what corresponds 
to the upper curve) appears as more restrictive than hav- 
ing one zero (which definitively is what corresponds to 
the lower curve). The generic four nonzero-probability 
case corresponds to g = 1, and the generic three nonzero- 
probability case corresponds to g = 0. Consequently, 
it appears as reasonable that the generic two nonzero- 
probability case would correspond to q — > — oo. This 
kind of paradoxal situation appears to need further dis- 
cussion to be clarified. The situation would probably be 
clear-cut if we had — but we do not! — , for arbitrary 
values of N arbitrary subsystems, the general answer to 
be associated with arbitrary q. This point thus remains 
open. 



a\ 


1 


2 




1 


p'^ + K 


p(l -p) - K 


P 


2 


P(l -P) - « 


{1-pf+K 


l-p 




P 


l-p 


1 



TABLE I: Joint probabilities for two binary subsystems A 
and B. The marginal probabilities are indicated as well. The 
correlation k and the probability p are such that all terms 
p^ + K,, p (1 — p) — K, and {l — p)^ + K remain within the interval 
[0, 1]. The case of independency corresponds to k = 0. 

It is worthy stressing a logical consequence of what has 
been said up to now. Unless g = 1, we cannot simultane- 
ously have equal probabilities for the set of joint probabil- 
ities of the allowed states and for the associated marginal 
probabilities. For example, for the q = case of two 
equal binary subsystems, equal probabilities for the joint 
set means = p^+^ = p^+^ = 1/3 and p4+^'= 0, 

hence p = 2/3 (which differs from 1 — p = 1/3), whereas 
equal probabilities for the marginal sets means p = 1/2, 



0.25- 
020- 
0.15- 
0.10- 
0.06- 
0.00- 



-0.0^ 
-0.10- 
-0.1^ 
-0.20- 



-0.25- 



/ a=1 


s^:— >+00 






\v q=0.5 yf 









FIG. 2: Curves k(p) which, for typical values of g, imply 
additivity of Sq. For — l/4<«;<0we have ^—k < p < 
1 - ^/ ^. Fo r < «: < 1/4 we have (1 - ^\ - 4k)/2 < p < 
(1 + VI - 4k)/2 . 



1/2. 



The q = 

~ /'21 



1 



hence +^ = ^ ^^3+^ - 
case simultaneously allows for p^j^^ 
P22^^ ==1/4 and p = 1 - p = 1/2. It is clear that, if 
such behavior (i.e., generic impossibility of equal proba- 
bilities for both joint and marginal ones simultaneously) 
persists up to the thermodynamic limit 00, it will 

have heavy consequences for the macroscopic statistics of 
the system. 



Let us recall at this stage that, depending on whether 
the subsystems that we are composing are or are not 
independent, it is Sbg or a different entropy which is ad- 
ditive. If they are correlated in the special form that has 
been illustrated above, and that persists up to the ther- 
modynamic limit iV — > 00, it is S'q with a specific value 
of q 7^ 1 which becomes extensive [l9l |. The correspond- 
ing statistical mechanics should consistently be based on 
Sq (or on a directly related one, such as say 5*2-^) and 
not anymore on Sbg- We may summarize this scenario 
through the following statement: Unless the composition 
law of the .subsystems is specified, the question whether an 
entropy (or some .similar quantity) is or is not extensive 
has no sense. Allow us a quick digression. The situa- 
tion is totally analogous to the quick or slow motion of a 
body. Ancient greeks considered motion to be an abso- 
lute property. It was not until Galileo that it was clearly 
perceived that motion has no sense unless the referential 
is specified. 
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3. The hint of the q- generalization of the Pascal triangle 

Very recently, Suyari and Tsukuda |23| used the con- 
cept of g'-product to consistently generalize various re- 
sults that are widely known in the context of BG sta- 
tistical mechanics and its mathematical structure. The 
generalization of the n\ (n-factorial) operation, as well as 
of the binomial (and even multinomial) coefficients was 
performed. As a corollary, the Pascal triangle itself was 
q-generalized as well. It is known that the coefficients 
of the n-th line of the Pascal triangle yield (after appro- 
priate centralization and rescaling) the Gaussian distri- 
bution in the n — > c» limit. This is well known to be 
a simple consequence of the standard central limit theo- 
rem as applied to independent binary variables. Suyari 
claims that, in the n ^ oo limit of the g-generalized Pas- 
cal triangle, what is obtained is precisely the g-Gaussian 
distributions! 



4- The hint of the scale- free networks 

The mathematical study of random networks (or ran- 
dom graphs) started many decades ago. But quite re- 
cently, it has acquired great interest due to the fact that 
such structures appear to be ubiquitous in physi cal, so- 
cial, internet and other complex phenomena (see |2lll2^ 
and references therein) . A central quantity of such struc- 
tures is the so-called connectivity or degree distribution, 
defined as the probability distribution of the number of 
links that are connected to the same site (or node); being 
more explicit, what one counts is the percentage of nodes 
that have a given number of links. For the important 
class of networks that are referred to as scale- free ones, 
this distribution is systematically found to be precisely a 
q-exponential. Many examples do exist in the literature. 
Three recent illustrations can be found in [2^ 123 . l25l |. 

These scale-invariant structures exhibit hubs and sub- 
hubs, and so on. If embedded into a d-dimensional space, 
they tend to have zero Lebesgue measure. They appear 
to be like fractals, which implies that most of the loca- 
tions of the d-dimensional space are forbidden (like the 
flights of any air company are expected to start and end 
only at the airports where that company operates, and 
not in any place of the territory!). 



5. The nonlinear dynamical hint 

We review at this point a few nonlinear dynamical re- 
sults which, although not directly related to the possible 
third central limit theorem we are seeking, provide what 
we consider important connections within the mathemat- 
ical structure which sustains noncxtcnsive statistical me- 
chanics. 

We focus on classical nonlinear dynamical systems, ei- 
ther conservative (e.g., Hamiltonian systems) or dissipa- 
tive. The basic ideas are quite general, but their presen- 



tation becomes easier if we illustrate them on a simple, 
one-dimensional, system. Let us address unimodal one- 
dimensional maps, such as for example the z-logistic map 
defined through 



xt+i = 1 - a\xt\'' 

■ ■ ; 2 > 1; < a < 2; -1 < xt < 1)(79) 

The sensitivity ^ to the initial conditions is defined as 
follows: 



(i = 0,l,2 



Ax{t) 



(80) 



where Ax(t) is the discrepancy at time t of two values of 
xq initially separated by Aa;(0). The typical behavior of 
f is given by 



(81) 



where Ai is referred to as the Lyapunov exponent (the 
subindex 1 will become clear in a few lines). The follow- 
ing ordinary differential equation is satisfied: d£,/dt = 

For most values of the control parameter a (see Eq. 
(79)), Ai is nonzero. When it is positive, we shall say 
that there is strong chaos. When it is negative, we have 
regular orbits as attractors (e.g., fixed points, cycles-2, 
cycles-3, cycles-4, etc). But there is an infinity of val- 
ues of a for which Ai vanishes. Examples are (i) all the 
values of a for which there is a doubling-period bifurca- 
tion from one cycle to its double; (ii) all the values of a 
for which there is a tangent bifurcation; (iii) all the val- 
ues of a for which there is an edge to chaos. All these 
cases are interesting, but by far the richest one is case 
(iii), referred to as weak chaos. Indeed, it is there that 
we expect the appearence of complexity (in physics, bi- 
ology, economics, linguistics, and elsewhere), since it is 
the frontier between considerable order (regular orbits) 
and considerable disorder (strong chaos). So, the follow- 
ing question arises: what is the behavior of ^{t) when Ai 
vanishes? The typical behavior is as follows: 



(82) 



where we have q > 1 (with Xq > 0) for the above cases (i) 
and (ii), and g < 1 (with A, > 0) for the above case (iii). 
The following ordinary differential equation is satisfied: 
d^/dt = Xq ^'J. The behavior (82) was conjectured in 1997 
[2g | and was proved in 2003 [23. For example, for z — 2, 
the first edge of chaos appears at a = Oc = 1.4011.... 
Its associated values are Qsen = 0.2445... and Ag,^„ — 
1/(1 — Qsen) — 1.3236..., where the subscript sen stands 
for sensitivity. 

Let us now address another very important quantity, 
namely the entropy production K per unit time. It is 
basically defined as follows. Take the admissible phase 
space (the interval [—1,1] for the variable x in the map 
(79), for instance), and make a partition of it in W little 
parts, identified through i = 1,2,...,M^. Then consider 
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Nq initial values of xq within one of the W intervals. As 
a function of time, the points within the initial window 
will possibly spread around, in such a way that we have 
N.i{t) points in the i-th interval iJ2Zi = ^o)- We 

define next the set of probabilities pi = Ni{t)/No (Vi), 
and finally we calculate the entropy Sq{t) through Eq. 
(9) for any chosen value of q. By construction we have 
Sq{0) — 0. We then define the entropy production per 
unit time as follows 



Kq — lim lim lim 

t—^OC W—^CXD Nq^OO 



t 



(83) 



The rate Ki coincides, in most cases, with the so called 
Kolmogorov-Sinai entropy rate. The rate Kq represents 
its ^-generalization. Consistently with the Pesin theo- 
rem, we expect, whenever Ai > 0, 



(84) 



This is indeed verified (see, for instance, |23|). For exam- 
ple, for (z, a) = (2, 2) in (79), we obtain Ki = Xi ^ In 2. 
What happens however at the edge of chaos (say at 
ac{z))7 Can we state something more informative than 
just Ki = Ai = 0? Yes, we can. It can be numerically 
shown [23 and analytically proved [s^ that 

^9.e„ = Ag,,„ . (85) 

This interesting relation (which refers in fact to upper 
bounds: see details in 1^3) basically generalizes the cele- 
brated Pesin theorem, and was also conjectured in 1997 
[26|. To be more specific, what happens is that there is a 
value of q (and only one), noted Qsen, such that Kq — 
for q > qsen, and Kq^ ooior q< qsen- And for q = qsen 



we obtain a finite value Ka 



with Ao 



When Ai > 0, qs 



which precisely coincides 
= 1 (strong chaos); when 



Ai = and the system is at the edge of chaos, qsen < 1 
and generically A,^^„ > (weak chaos). For example, 
for the edge of chaos of the universality class to which 
the map (79) belongs, qsen increases from —00 to slightly 
below unity when z increases from 1 to cx). 

What happens when our nonlinear dynamical system 
has more than one, say /i nonnegative Lyapunov expo- 
nents {A^^™-*} (to = 1,2, ...,/i)? From the Pesin theorem, 
we certainly expect 



^1 = E ^i' 



(m) 



(86) 



What happens then if the system is at an edge of chaos, 
where A^™"* = , Vto? We expect js^ the following (con- 
jectural) behavior for the Lebesgue measure ^ associated 
with the dynamically expanding directions: 



(87) 



where the (~) sign might become just {—) under some 
simplifying hypothesis (like orthogonality of the direc- 
tions along which the expansions associated with the 
X\ s occur). 



Since we essentially expect, for t — > cx). 



(the subscript e stands for entropy), we possibly have the 
following relation: 



So. 



In, 



r ' 

[n 



m— 1 



(89) 



Using the definition (83) we obtain the following inter- 
esting relation: 



ln„ 



Kr, 



lim ■ 

t — >oo 



t 



(90) 



Two important cases must be distinguished, namely 
strong and weak chaos. Strong chaos corresponds to 

A^™^ > , Vm. In this case, we have qe = qiTn = 1 , Vm, 
hence Eq. (90) straightforwardly recovers relation (86), 
consistently with the Pesin theorem. Weak chaos corre- 
sponds to A^"*^ = , Vto. In this case, by focusing on the 
t ^ 00 asymptotic region, we have 



K„ 



m=l - qsen )^ J 



1 - 

,(i-'je)i::;=i 



X lim 

t — ^oo 



t 



(91) 



consequently (since qe must be chosen so that Kq^ is fi- 
nite) 



and 



Kq,= 



1 - ge 



11™= 1 



1 

E -, _ (m) ' 
m=l qsen 



(92) 



(m)^^ , (m) 



[ 1 ~ qsen ) A (,„) 



1 - ge 



(93) 



Eq. (93) can be rewritcn in a more symmetric form, 
namely 



[{l-qe)KqJ^=l[\{l-qiT^)X% 

L Qsen 



(94) 



If ^ — 1, Eq. (92) recovers qe = qsen, and Eq. (94) 
recovers Eq. (85). 

Another interesting particular case is when qi^J — 
qsen (Vm), and A*"™', — Xq^^^ (Vto). From Eqs. (92) and 
(94) it follows then 



1 - ge 



1 - 



(95) 
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and 



(96) 



An example which verifies Eq. (95) can be found in |3lj |. 

One more particular case which is interesting is when 
the expansion is linear in time for all directions, i.e., 
(Vm). We then have, from Eqs. (92) and (94), 



1 

1 



and 



r ^ 

[n 



1/m 



(97) 



(98) 



Relation (97) has already emerged in the literature 
through various forms. A first example concerns /i = 1, 
hence we expect qe = 0. This is precisely what can be ver- 
ified for the Casati-Prosen triangle map 33], which 
is a two-dimensional, conservative (hence /i = 1), mix- 
ing, ergodic one, with linear instability. In addition to 
Qe — 0, it has been verified that K^. = Aq, in accordance 
with Eq. (98). A second example is the ^-dimensional 
lattice Lotka-Volterra, for which it has precisely been 
verified [3 IHl the result (97). A third example is a 
specific c?-dimensional Boltzmann latice model 36]. Its 
Hamiltonian-like behavior leads to /i = d/2, which, from 
Eq. (97), implies Op = l~2/d. It is precisely this result 
that the authors ISo] have obtained by imposing Galilean 
invariance to the dynamical equations of their model. 

It is clear that the interesting (and possibly exact) re- 
lations (92) and (94) remain to be proved. At the present 
stage they constitute but conjectures. 



6. The empirical hint 

Last but not least, let us present a very pragmatic, 
epistemological-like, reason. There are, in the literature, 
already so many natural and artificial systems (and their 
number constantly increases) whose central quantities are 
well fitted by q-exponentials and q-Gaussians, that one 
is compelled to believe that only a limit theorem could 
explain such an ubiquity. 



V. CONCLUSIONS 

In this final section, let us summarize the scenario 
within which we are working. It appears that the entropic 
form S{N; {pi}) to be used for constructing a statistical 
mechanics that is naturally compatible with usual macro- 
scopic thermodynamics should be generically extensive, 
i.e., such that 



< lim — < oo . 



(99) 



the equality occuring only for the case of certainty or 
(thermodynamically) close to it (the typical example is 
a system in themal equilibrium at zero temperature). 
Within this (axiomatic) viewpoint, we can distinguish 
the following cases: 

(i) The N subsystems (typically physical elements such 
as particles free to move translationally and/or rotation- 
ally, localized spins, and similar entities) are (probabilis- 
tically) independent (the paradigmatic case is that of 
ideal gases; their total energy is strictly proportional to 
N, i.e., additive, hence extensive). Then we must use 
Sbg since SBciN; {pt}) cx N (Vfc}). 

(ii) The N subsystems are locally correlated (the 
paradigmatic case is that of s/iort-range-interacting 
many-body Hamiltonian systems; their total energy is 
asymptotically proportional to N , i.e., it is extensive once 
again). Then, once again, wc must use Sbg since it sat- 
ifiesEq. (99), Vfe}. 

(iii) The N subsystems are globally correlated in the 
special manner addressed in this paper (the paradigmatic 
case appears to be that of lora(7-range-interacting many- 
body Hamiltonian systems; their total energy asymptot- 
ically increases with N faster than linearly, i.e., it is 
nonextensive) . Then, we must use Sq for that unique 
value of q which guarantees Eq. (99), V{pi}. 

(iv) The N subsystems are globally correlated in a man- 
ner which is more complex (or just as complex but in a 
different manner) than the one addressed in this paper. 
Then we must use an entropy which is not included in the 
family Sq for any value of q. Such an entropy would have 
to be either more general than Sq (see for instance for 
the so called superstatistics '38]), or just of a completely 
different type (see for instance 44] ) . 

The cases (i) and (ii) , which we may call simple in the 
sense of plectics j43 ]. belong to the world within which 
the concepts used in BG statistical mechanics have, since 
more than one century, been profusely shown to be the 
appropriate ones. 

The cases (iii) and (iv), which we may call complex in 
the sense of plectics jij , belong to the world within which 
the concepts used in nonextensive statistical mechanics 
(as well as in its possible generalizations or alternatives) 
have been shown and keep being shown (since more than 
one decade, by now) to be the appropriate ones. 

The cause for a system to be simple or complex in the 
above sense lies basically (see 0,|43) on its microscopic 
dynamics in the full space of microscopic possibilities 
(Gibbs' F-space for many-body Hamiltonian systems). It 
is believed to be so because it is this dynamics which is 
expected to determine the possible persistent correlations 
that would be responsible for the geometrical structure 
within which the system tends to "live" , given its initial 
conditions. 

If its elementary nonlinear dynamics is controlled by 
strong sensitivity to the initial conditions (i.e., at least 
one positive Lyapunov exponent if the system is a clas- 
sical one), we expect the system to be simple. For vir- 
tually any initial condition, the system quickly visits the 
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neighborhood of virtually all possible states; it does so in 
such a way that the probability of any small part of the 
(macroscopically) admissible full space asymptotically 
becomes proportional to its size (Lebesgue measure) . Its 
main time-dependent functions (typically) exponentially 
depend on time, the number of allowed stated exponen- 
tially increases with N, and, if the system is Haniiltonian, 
its stationary state (called thermal equilibrium) is charac- 
terized by an energy distribution given by the celebrated 
BG weight. Summarizing, its paradigmatic ordinary dif- 
ferential equation is dy/dx cx y. 

If the elementary nonlinear dynamics of the system is 
nonregular and controlled by weak sensitivity to the ini- 
tial conditions (i.e., no positive Lyapunov exponents if 
the system is a classical one), we expect the system to 
be complex. For given initial conditions, the system es- 
sentially visits a network of states (a scale-free network 
for many if not all g-systems) whose typical Lebesgue 
measure is zero. The particular network depends from 
the initial conditions and is highly inhomogeneous (like, 
as mentioned before, the network of airports on which a 
specific air company operates), but its geometry (both 
topology and metrics) is basically the same for virtu- 
ally all initial conditions. To recover a homogeneous oc- 
cupancy of the full space we are obliged to make av- 
erages over all the initial conditions (whereas no such 
thing is necessary for simple systems). The main time- 
dependent functions of the system typically depend on 
time slower than exponentially (typically like power-laws, 
more precisely like q-exponentials for most if not all the 
q-systems), the number of allowed states increases with 



N like a power-law, and, if the system is Hamiltonian, its 
stationary or quasi-stationary (metastable) state (out of 
thermal equilibrium) is expected to be characterized by 
an energy distribution given by a g-exponential weight, 
which asymptotically approaches a power-law for high 
energies. Summarizing, its paradigmatic ordinary differ- 
ential equation is dy/dx cx y'' . 

The standard central limit theorem plays a crucial role 
for the simple systems. We expect a similar theorem to 
exist for the complex systems of the g-class. The proof of 
such a theorem would be priceless. 

In addition to the above, we have presented here more 
two conjectures (Eqs. (92) and (94)) concerning the en- 
tropy production per unit time for a nonlinear dynamical 
system at the edge of chaos, having n vanishing Lyapunov 
exponents. The entropy Sq which increases linearly with 
time (when the system is exploring its phase space) ap- 
pears to be that whose entropic index is Qe as given by 
Eq. (92) [sj- The associated entropy production i^Tg^ per 
unit time appears to be as given by Eq. (94). The proof 
of these two connected conjectures (including, naturally, 
the precise conditions for their validity) also remains as 
a mathematically open problem. 

Acknowledgments 

We acknowledge useful discussions with S. Abe, C. 
Anteneodo, F. Baldovin, S. Curilef, M. Gell-Mann, R. 
Hersh, L. Moyano, C. Pagani and Y. Sato, as well as 
partial financial support by Pronex/MCT, Faperj and 
CNPq (Brazilian agencies). 



[1] L. Boltzmann, Weitere Studien uber das 
Wdrm,egleichgewicht unter Gas molekulen [Further Stud- 
ies on Thermal Equilibrium Between Gas Molecules], 
Wien, Ber. 66, 275 (1872), 

[2] J.W. Gibbs, Elementary Principles m Statistical Me- 
chanics - Developed with Especial Reference to the Ratio- 
nal Foundation of Thermodynamics (C. Scribner's Sons, 
New York, 1902; Yale University Press, New Haven, 1948; 
OX Bow Press, Woodbridge, Connecticut, 1981). In page 
35: "In treating of the canonical distribution, we shall al- 
ways suppose the multiple integral in equation (92) [the 
partition function, as we call it nowadays] to have a finite 
value, as otherwise the coefficient of probability vanishes, 
and the law of distribution becomes illusory. This will ex- 
clude certain cases, but not such apparently, as will affect 
the value of our results with respect to their bearing on 
thermodynamics. It will exclude, for instance, cases in 
which the system or parts of it can be distributed in un- 
Umited space [...]. It also excludes many cases in which 
the energy can decrease without limit, as when the sys- 
tem contains material points which attract one another 
inversely as the squares of their distances. [...]. For the 
purposes of a general discussion, it is sufficient to call 
attention to the assumption implicitly involved in the 
formula (92)." 

[3] A. Einstein, Theorie der Opaleszenz von homogenen 



Flussigkeiten und Flussigkeitsgemischen in der Ndhe des 
kritischen Zustandes, Annalen der Physik 33, 1275 
(1910). ["Usually W is put equal to the number of com- 
plexions... In order to calculate W, one needs a com- 
plete (molecular-mechanical) theory of the system un- 
der consideration. Therefore it is dubious whether the 
Boltzmann principle has any meaning without a complete 
molecular-mechanical theory or some other theory which 
describes the elementary processes. S = log W-\- const, 
seems without content, from a phenomenological point of 
view, without giving in addition such an Elementarthe- 
orie." (Translation: Abraham Pais, Subtle is the Lord..., 
Oxford University Press, 1982)]. 

[4] M. Gell-Mann and C. Tsallis, eds., Nonextensive En- 
tropy - Interdisciplinary Applications (Oxford University 
Press, New York, 2004); H.L. Swinney and C. Tsallis, 
eds.. Anomalous Distributions, Nonlinear Dynamics and 
Nonextensivity, Physica D 193 (2004). 

[5] C. Tsallis, Possible generalization of Boltzmann- Gibbs 
statistics, J. Stat. Phys. 52, 479 (1988); for full bibli- 
ography see http://tsallis.cat.cbpf.br/biblio.htm 

[6] E.M.F. Curado and C. Tsallis, Generalized statistical me- 
chanics: connection with thermodynamics, J. Phys. A 24, 
L69 (1991); Corrigenda: 24, 3187 (1991) and 25, 1019 
(1992). 

[7] C. Tsallis, R.S. Mendes and A.R. Plastino, The role 



13 



[8] 



[9] 
[10] 

[11] 
[12] 

[13] 
[14] 

[15] 



[16] 

[17] 
[18] 



of constraints within generalized nonextensive statistics, 
Physica A 261, 534 (1998). 

F. Takens, in Structures in Dynamics - Finite Dimen- 
sional Deterministic Studies, eds. H.W. Broer, F. Du- 
mortier, S.J. van Strien and F. Takens, (North-Holland, 
Amsterdam, 1991). In page 253: "The values of Pi are 
determined by the following dogma : if the energy of 
the system in the i*'' state is Ei and if the tempera- 
ture of the system is T then: pi = e~^^^^'^ /Z{T), where 
Z{T) = e-^"'*'^/Z(r), (this last constant is taken so 
that YliPi ~ 1)- This choice of pi is called the Gibbs dis- 
tribution. We shall give no justification for this dogma; 
even a physicist like Ruelle disposes of this question as 
"deep and incompletely clarified"." 

C. Tsallis, What are the numbers that experiments pro- 
vide?, Quimica Nova 17, 468 (1994). 

D. Prato and C. Tsallis, Nonextensive foundation of Levy 
distributions, Phys. Rev. E 60, 2398 (1999). 

A.M.C. de Souza and C. Tsallis, Student's t- and r- 
distributions: Unified derivation from an entropic vari- 
ational principle, Physica A 236, 52 (1997). 
A. Einstein, Uber die von der molekular kinetischen The- 
orie der Warme gefordete Bewegung von in ruhenden 
Flussigkeiten suspendierten Teilchen, Annalen der Physik 
17, 549 (1905). 

A.R. Plastino and A. Plastino, Non-extensive statistical 
mechanics and generalized Fokker-Planck equation, Phys- 
ica A 222, 347 (1995)1. 

C. Tsallis and D.J. Bukman, Anomalous diffusion in the 
presence of external forces: exact time- dependent solu- 
tions and their thermostatistical basis, Phys. Rev. E 54, 
R2197 (1996). 

E. M.F. Curado and F.D. Nobre, Derivation of nonlin- 
ear Fokker-Planck equations by means of approximations 
to the master equation, Phys. Rev. E 67, 021107 (2003); 

F. D. Nobre, E.M.F. Curado and G. Rowlands, A proce- 
dure for obtaining general nonlinear Fokker-Planck equa- 
tions, Physica A 334, 109 (2004). 

L. Nivanen, A. Le Mehaute and Q.A. Wang, General- 
ized algebra within a nonextensive statistics. Rep. Math. 
Phys. 52, 437 (2003); E.P. Borges, A possible deformed 
algebra and calculus inspired in nonextensive thermo- 
statistics, Physica A 340, 95 (2004). 
A simple generalization also covering possibly negative 



[19] 



real numbers X and Y is the following one: X < 



Y 



[sign{X)][sign{Y)][\X\ 



+ y 



- 1 



1/(1-9) 



in Complexity, Metastbility and Nonextensivity, eds. C. 
Beck, A. Rapisarda and C. Tsallis (World Scientific, Sin- 
gapore, 2005), to appear cond-mat/0409631 ; Y. Sato 
and C. Tsallis, On the extensivity of the entropy Sq for 
systems with N < 3 specially correlated subsystems. In- 
ternational Journal of Bifurcation and Chaos, to appear: 
Proceedings of the Summer School and Conference on 
Complexity in Science and Society, ed. T. Bountis (Pa- 
tras and Olympia, 14-26 July, 2004) cond-mat/0411073 . 
To avoid ambiguities let us recall that an entropy S(N) 
of a system composed of A'^ subsystems is said additive 
if, under a given composition law for the N subsystems, 
S{N) oc A^. It is said extensive if limjv^cx) S{N)/N is fi- 
nite. Consequently, if an entropy is additive, it also is 
extensive. The opposite is not true. The usual thermody- 
namic requirement for entropy is to be extensive (i.e., 



[20] 



[21] 
[22] 
[23] 



[24] 

[25] 
[26] 
[27] 



[28] 
[29] 
[30] 

[31] 



asymptotically additive), not necessarily additive (i.e., 
strictly additive). So, if the subsystems are independent, 
we must use Sbg since it is additive, hence extensive. If 
the subsystems are locally correlated (the paradigmatic 
case is that of s/iori-range-interacting many-body Hamil- 
tonian systems), we must also use Sbg since, although 
not being strictly additive, it is extensive. If the sub- 
systems are globally correlated (the paradigmatic case 
appears to be that of tong^range-interacting many-body 
Hamiltonian systems), then we must use another entropy 
as the basis for constructing a statistical mechanics. If 
this global correlation is of the type focused on in this 
paper, then Sq with the appropriate value of q is exten- 
sive, and is the one to be used. If the global correlation 
is of some different nature, even more complex than the 
one focused on in this paper, a still different entropic 
form might be necessary. 

H. Suyari and M. Tsukada, Law of error in Tsallis statis- 
tics, IEEE Transaction on Information Theory (2004), 
in press cond- mat/0401540 ; H. Suyari, q-Stirlmg's for- 
mula in Tsallis statistics, cond-mat/0401541 H. Suyari, 
Mathematical structure derived from the q-multmomial 
coefficient in Tsallis statistics, cond-mat/0401546 
D.J. Watts and S.H. Strogatz, Collective dynamics of 
small-world networks, Nature 393, 440 (1998). 
R. Albert and A.-L. Barabasi, Statistical mechanics of 
complex networks, Rev. Mod. Phys. 74, 47 (2002). 
R. Albert and A.-L. Barabasi, Topology of evolving net- 
works: Local events and universality, Phys. Rev. Lett. 85, 
5234 (2000); the analytically exact distribution which is 
found in this paper is, although writen in a slightly differ- 
ent form, a q-exponential with q — , where 

' ^ ^ ^ m(3 — 2r) + l— p— r ' 

{m,p,r) are parameters of the specific growth model. 
J.P.K. Doye, Network topology of a potential energy land- 
scape: A static scale-free network,Phys. Rev. Lett. 88, 
238701 (2002); the distributions that are presented in 
this paper can be shown to be numerically very close to 
a q-exponential with q ~ 2. 

D.J.B. Soares, C. Tsallis, A.M. Mariz and L.R. da Silva, 
Preferential attachment growth model and nonextensive 
statistical mechanics, cond- mat/0410459 
C. Tsallis, A.R. Plastino and W.-M. Zheng, Power-law 
sensitivity to initial conditions - New entropic represen- 
tation. Chaos, Solitons and Fractals 8, 885 (1997). 
F. Baldovin and A. Robledo, Sensitivity to initial condi- 
tions at bifurcations in one- dimensional nonlinear maps: 
Rigorous nonextensive solutions, Europhys. Lett. 60, 
518 (2002); F. Baldovin and A. Robledo, Universal 
renormalization- group dynamics at the onset of chaos 
in logistic maps and nonextensive statistical mechanics, 
Phys. Rev. E 66, R045104 (2002). 

V. Latora and M. Baranger, Kolmogorov-Sinai entropy 
rate versus physical entropy, Phys. Rev. Lett. 82, 520 
(1999). 

V. Latora, M. Baranger, A. Rapisarda and C. Tsallis, 
The rate of entropy increase at the edge of chaos, Phys. 
Lett. A 273, 97 (2000). 

F. Baldovin and A. Robledo, Nonextensive Pesm iden- 
tity. Exact renormalization group analytical results for 
the dynamics at the edge of chaos of the logistic map., 
Phys. Rev. E 69, 045202(R) (2004). 

G. F.J. Ananos, F. Baldovin and C. Tsallis, Anoma- 
lous sensitivity to initial conditions and entropy pro- 



14 



duction in standard maps: Nonextensive approach, 

cond-mat/0403656 
[32] G. Casati, C. Tsallis and F. Baldovin, Linear instability 

and statistical laws of physics, in preparation. 
[33] G. Casati and T. Prosen, Mixing property of triangular [39] 

billiards, Phys. Rev. Lett. 83, 4729 (1999); G. Casati and 

T. Prosen, Triangle map: A model of quantum chaos, 

Phys. Rev. Lett. 85, 4261 (2000). 
[34] G.A. Tsekouras, A. Provata and C. Tsallis, Nonextensiv- 

ity of the cyclic lattice Lotka Volterra model, Phys. Rev. 

E 69, 016120 (2004). [40] 
[35] C. Anteneodo, Entropy production in the cyclic lattice 

Lotka- Volterra model, Eur. Phys. J. B (2004), in press 

cond-mat/0402248 . [41] 
[36] B.M. Boghosian, P.J. Love, P.V. Coveney, I.V. Karlin, S. 

Succi and J. Yepez, Galilean-invariant lattice-Boltzmann 

models with H-theorem, Phys. Rev. E 68, 025103(R) 

(2003). 

[37] C. Tsallis and A.M.C. Souza, Constructing a statistical 

mechanics for Beck-Cohen superstatistics, Phys. Rev. E [42] 
67, 026106 (2003); see also A.M.C. Souza and C. Tsallis, 
Stability of the entropy for superstatistics, Phys. Lett. A 
319, 273 (2003); A.M.C. Souza and C. Tsallis, Stability 
analysis of the entropy for superstatistics, Physica A 342, 
132 (2004). [43] 

[38] C. Beck and E.G.D. Cohen, Superstatistics, Physica A 
322, 267 (2003); see also C. Beck, Superstatistics: The- 
ory and Applications, in Nonadditive entropy and nonex- 
tensive statistical mechanics, ed. M. Sugiyama, Contin- 
uum Mechanics and Thermodynamics 16, 293 (Springer- 



Verlag, Heidelberg, 2004); E.G.D. Cohen, Superstatis- 
tics, in Anomalous Distributions, Nonlinear Dynamics 
and Nonextensivity, eds. H.L. Swinney and C. Tsallis, 
Physica D 193, 35 (2004). 

E.M.F. Curado, General aspects of the thermodynami- 
cal formalism, in Nonextensive Statistical Mechanics and 
Thermodynamics, eds. S.R.A. Salinas and C. Tsallis, 
Braz. J. Phys. 29, 36 (1999); E.M.F. Curado and F.D. 
Nobre, On the stability of analytic entropic forms, Phys- 
ica A 335, 94 (2004). 

C. Anteneodo and A.R. Plastino, Maximum entropy ap- 
proach to stretched exponential probability distributions, 
J. Phys. A 32, 1089 (1999). 

G. Kaniadakis, Nonlinear kinetics underlying generalized 
statistics, Physica A 296, 405 (2001); G. Kaniadakis, 
Statistical mechanics in the context of special relativity, 
Phys. Rev. E 66, 0561125 (2002); G. Kaniadakis, M. 
Lissia and A.M. Scarfone, Deformed logarithms and en- 
tropies, Physica A 340, 41 (2004). 

M. Gell-Mann, The Quark and the Jaguar - Adventures 
in the Simple and the Complex, (W.H. Freeman, New 
York, 1994). Plectics is the denomination suggested by 
M. Gell-Mann for the emerging science of the so called 
complex systems. 

E.G.D. Cohen, Boltzmann and Einstein: Statistics and 
dynamics - An unsolved problem, Boltzmann Award Lec- 
ture at STATPHYS 2004, 22nd lUPAP International 
Conference on Statistical Physics (Bangalore, 4 to 9 July 
2004), Pramana (2004), in press. 



